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Executive  Summaiy 


The  objective  of  this  program  was  the  development  of  a  truly  global  near  real-time 
ionospheric  specification  model  (PRISM;  Parameterized  Real-time  Ionospheric  Specification 
Model).  PRISM  is  based  on  parameterized  physical  models  rather  than  statistical  or 
climatological  models.  It  uses  near  real  time  data  from  ground  based  sources  (DISS 
bottomside  soundings)  and  satellite  bases  sources  (DMSP;  SSIES  in  situ  plasma 
measurements,  SSJ/4  precipitating  ion  and  electron  measurements,  and  electron  density 
profiles  derived  from  SSUSI  and  SSULI  images)  to  adjust  the  parameterized  models. 
Previous  specification  models  were  based  on  a  small  number  of  "global"  parameters  (e.g., 
effective  sunspot  number)  determined  from  a  limited  data  set.  The  adjustment  parameters 
in  PRISM  vary  with  location  and  can  be  determined  from  a  wide  variety  of  data 
representing  most  of  the  new  near  real-time  data  that  will  become  available  to  the  Space 
Forecast  Center  during  this  decade. 

Four  separate  physical  models  are  used  in  the  development  of  PRISM:  a  low  latitude 
F  layer  model,  a  midlatitude  F  layer  model,  a  combined  low  and  midlatitude  E  layer 
model,  and  a  high  latitude  E  and  F  layer  model.  To  achieve  sufficient  computational 
speed,  the  models  were  parameterized  in  terms  of  geophysical  parameters.  The 

parameterization  process  involved  the  production  of  "databases"  for  various  values  of  the 
geophysical  parameters  and  the  generation  of  semi-analytical  representations  of  the 
databases.  The  algorithm  for  adjusting  the  parameterized  models  on  the  basis  of  real  time 
data  is  described  in  this  report.  The  source  code,  an  input  data  set,  and  the  associated 
output  data  stream,  all  in  machine  readable  form,  were  delivered  to  the  Air  Weather 
Service  at  the  same  time  as  this  report. 

Some  data  has  been  acquired  and  more  data  is  being  acquired  for  testing  and 

validating  PRISM.  Actual  validation  will  begin  immediately  so  that  it  can  be  completed 

before  PRISM  is  released  for  transition  to  operational  status.  A 


IV 


supplemental  report  describing  the  testing  and  validation  process,  PRISM  performance,  and 
a  complete  error  propagation  analysis  of  PRISM,  will  be  delivered  at  the  same  time. 

Anticipated  future  enhancements  to  PRISM  include  the  use  of  TISS  data  and  the 
inclusion  of  H'*’  to  extend  PRISM  to  the  plasmapause. 


V 


Section  1:  Introduction 


This  report  describes  Version  1.0  of  PRISM  as  delivered  to  the  Air  Weather  Serv'ice. 
This  section  describes  the  development  objectives  and  our  approach  to  achieving  them. 
Section  3  describes  the  physical  models  on  which  PRISM  is  based.  Section  4  describes 
how  the  physical  models  were  parameterized.  Section  5  describes  how  the  real  time  data 
are  u.sed  to  adjust  the  parameterized  model.  Section  6  describes  the  structure  of  the  code, 
its  input  and  output,  and  describes  the  results  of  some  test  runs  using  simulated  data. 
References  cited  in  the  text  are  listed  in  Section  7. 

Section  1. 1  Objectives 

The  primary  objective  of  this  effort  is  the  development  of  an  algorithm  for  using  near 
real  time  satellite  and  ground  based  data  to  provide  a  near  real  time  specification  of  the 
global  ionosphere.  The  data  that  are  to  be  used  include 

(1)  bottomside  soundings  from  the  DISS  network, 

(2)  Total  Electron  Content  (TEC)  data  from  the  TISS  network, 

(3)  in  situ  plasma  data  (densities,  temperatures,  and  drift  velocities)  from  the  SSIES 
instrument  on  DMSP  satellites, 

(4)  auroral  electron  and  ion  fluxes  from  the  SSJ/4  instrument  on  DMSP  satellites,  and 

(5)  electron  density  profile  information  deduced  from  the  ultraviolet  instruments 
(SSUSI  and  SSULI)  expected  to  be  flown  on  DMSP  satellites. 

Unfortunately,  complications  in  making  use  of  TISS  data  (described  below)  prevented  the 
use  of  TISS  data  in  Version  1.0.  We  anticipate  correcting  this  deficiency  in  the  future. 
The  other  four  data  categories  are  used  in  Version  1.0. 

The  need  for  a  global  specification  of  the  state  of  the  ionosphere  is  twofold.  First, 
there  are  operational  systems  that  need  to  correct  for  ionospheric  effects  in  real  time,  or 
that  have  operational  parameters  that  are  affected  by  the  ionosphere  and  must  be  adjusted 
in  real  time.  Further,  the  operation  of  these  systems  could  be  optimized  if  accurate 
forecasts  of  ionospheric  conditions  were  available  because  this  allows  the  operational 
parameters  to  be  chosen  ahead  of  time.  Since  any  ionospheric  forecast  algorithm  will 
require  an  accurate  specification  of  the  current  state  of  the  ionosphere  as  an  initial 
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condition,  this  is  the  second  reason  an  accurate  ionospheric  specification  algorithm  is 
urgently  needed. 

Section  1.2  Approach 

Ideally,  the  specification  of  the  current  state  of  the  ionosphere  would  be  obtained 
directly  from  real  time  observations  from  a  dense  network  of  satellite  and  ground  based 
instruments.  Unfortunately,  the  complexity  of  the  ionospheric  system  realistically  precludes 
the  deployment  of  a  sufficiently  dense  network  of  observing  instruments.  Therefore,  any 
ionospheric  specification  algorithm  must  be  based  on  an  ionospheric  model  with  parameters 
that  can  be  adjusted  on  the  basis  of  near  real  time  data.  Two  approaches  are  possible:  (1) 
statistical  or  climatological  models  and  (2)  numerical  simulations  based  on  physical  models. 

For  reasons  described  below,  we  have  chosen  the  second  approach  (physical  models). 
However,  practical  considerations  (primarily  computational  speed)  dictate  that  the 
algorithms  implemented  at  the  Space  Forecast  Center  be  based  on  parameterized  versions  of 
the  physical  models.  Thus,  we  call  the  global  algorithm  PRISM  for  Parameterized 
Real-time  Ionospheric  Specification  Model. 

Section  1.2.1  Motivation 

We  feel  strongly  that  a  comprehensive  physical  model  of  ionospheric  proces.ses  can 
produce  more  accurate  specifications  and  forecasts  than  can  statistical  or  climatological 
models.  The  causal  relationship  between  easily  monitored  solar  and  geophysical  parameters 
(e.g.,  K^,  etc.)  and  particular  ionospheric  configuration  is  very  complex.  Any 

organization  of  historical  ionospheric  data  inevitably  averages  over  a  variety  of 
configurations  corresponding  to  similar  values  of  the  chosen  .set  of  solar-geophysical 
parameters  (usually  only  one  or  two).  The  result  is  that  spatial  structure  tends  to  be 
smeared  out  or  smoothed  over  and,  therefore,  unrepresentative  of  the  instantaneous 
ionosphere.  If  a  physical  model  contains  all  of  the  relevant  physics,  then  it  will  produce 
more  realistic  representations  of  instantaneous  ionospheric  structure.  However,  there  is  a 
difference  between  a  realistic  representation  and  an  accurate  one. 


In  order  to  accurately  simulate  a  time  dependent  phenomenon  like  the  ionosphere,  a 
physical  model  needs  an  accurate  specification  of  the  initial  conditions  and  an  accurate 
representation  of  the  energy  and  momentum  flux  at  the  boundaries.  For  the  purpo.ses  of 
providing  a  specification  model,  it  is  the  energy  and  momentum  input  that  is  crucial.  If 
the  model  is  run  long  enough,  the  effects  of  the  initial  conditions  are  lost  and  the  present 
state  of  the  model  depends  only  on  the  recent  history  of  the  energy  and  momentum  input. 
These  include  the  solar  EUV  (the  primary  source  of  ionization  outside  the  auroral  zone), 
high  latitude  heating  of  the  thermosphere  (which  affects  the  global  circulation  of  the 
thermosphere),  high  latitude  convection,  and  low  latitude  dynamo  electric  fields.  While  the 
temporal  and  spatial  resolution  of  the  observations  of  these  quantities  are  expected  to 
improve  in  the  future,  they  will  probably  always  be  insufficient  to  allow  accurate 
ionospheric  simulation  of  the  ionosphere  without  additional  data.  As  a  practical  matter, 
ionospheric  simulations  must  be,  and  will  remain,  iterative  in  nature.  The  energy  and 
momentum  input  parameters  are  adjusted  until  the  simulation  agrees  with  observations  of 
ionospheric  parameters  to  some  level  of  accuracy. 

A  practical  consequence  of  this  situation  is  that  the  production  of  an  accurate 
ionospheric  specification  based  on  a  numerical  simulation  requires  that  the  physical  model 
be  run  several  times.  The  requirement  that  the  model  cover  a  sufficient  time  period  to 
allow  the  transient  effects  of  the  initial  conditions  to  damp  out  implies  that  the  physical 
model  must  run  much  faster  than  the  system  it  is  simulating.  At  the  present  time,  given 
practical  limits  on  the  available  computing  power,  this  is  not  possible.  Consequently,  we 
have  adopted  a  modified  approach  in  which  the  physical  models  are  parameterized  in  terms 
of  solar  and  geophysical  parameters.  It  is  these  parameterized  models  rather  than  the 
original  physical  models  which  are  to  be  adjusted  according  to  the  real  time  ionospheric 
data. 


There  is  a  superficial  similarity  between  our  approach  and  a  climatological  approach. 
The  difference,  however,  is  that  we  begin  with  a  more  realistic  representation  of  the  spatial 
structure  of  the  ionosphere  than  climatological  models  can  provide.  The  parameter 
adjustment  process  should  not  compromise  this  advantage.  In  the  future,  as  more  powerful 
computers  and  more  efficient  model  algorithms  become  available,  the  parameterized  models 
can  be  replaced  by  actual  physical  models  to  produce  more  accurate  specifications. 
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Section  1.2.2  The  Parameterization  Process 


The  physical  models  are  parameterized  in  two  steps.  First,  the  physical  models  are 
run  for  a  selected  set  of  values  of  the  solar  and  geophysical  parameters  (e.g.,  season,  solar 
activity,  and  magnetic  activity)  to  produce  a  set  of  "databases".  Each  database  consists  of 
ion  density  profiles  for  a  24  hour  period  over  the  range  of  latitudes  and  longitudes  for 
which  the  model  is  applicable.  Separate  models  are  used  for  the  low,  middle,  and  high 
latitudes  with  care  taken  to  ensure  smooth  transitions  at  region  boundaries.  Second, 
semi-analytical  representations  of  each  database  are  produced.  The  altitude  profiles  are 
represented  by  linear  combinations  of  Empirical  Orthonormal  Functions  (EOF’s,  see  Section 
3.2.1)  that  are  generated  for  each  database.  The  longitude  variation  of  the  coefficients  of 
the  EOFs  is  represented  by  a  truncated  Fourier  series,  and  the  latitude  variation  of  the 
Fourier  coefficients  is  represented  by  a  linear  combination  of  orthogonal  polynomials 
(Section  3.2.2).  Because  the  EOFs  are  tabulated  functions,  the  representation  is  termed 
"semi-analytic." 


Section  1.2.3  The  Real  Time  Adjustment  Algorithm 

The  semi-analytic  representation  of  each  ionospheric  "database"  allows  the  results  of 
the  physical  model  runs  to  be  reconstructed  quickly  and  accurately.  However,  each 
database  corresponds  to  a  specific  combination  of  energy  and  momentum  input  conditions, 
and  only  rarely  will  the  actual  conditions  coincide  with  any  one  of  the  databa.ses.  Because 
of  the  nonlinear  nature  of  the  ionosphere,  and  because  each  database  is  represented  by  a 
separate  set  of  basis  functions  (EOF's),  the  Version  1.0  of  PRISM  does  not  interpolate 
between  databases.  Instead,  the  database  corresponding  most  closely  to  the  actual 
conditions  is  adjusted  on  the  basis  of  the  available  near  real  time  data.  The  details  of  the 
adjustment  process  are  described  in  Section  4  but  will  be  summarized  here.  The  high 
latitude  region  is  handled  differently  than  the  low  and  midlatitude  regions. 

Low  and  Middle  Latitudes.  The  corrections  to  the  nominal  electron  density  profile  will  be 
characterized  by  eight  parameters  that  are  functions  of  latitude  and  longitude. 

(1)  critical  frequencies, 

(3)  Ah  Fr.  and  (4)  Ah  E,  corrections  to  the  heights  of  the  E  and  F  layer  peaks. 
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(5)  ^top'  topside  electron  density  and  scale  height 

at  the  DMSP  altitude 

(7)  A/Vj  and  (8)  AN^,  corrections  to  the  electron  density  at  8  km  and  16  km  below 
the  E  layer  peak,  respectively. 

These  parameters  represent  a  correction  to  the  parameterized  model  as  determined  by 
comparison  with  near  real  time  data.  Thus,  the  nominal  value  of  each  of  these  parameters 
is  zero,  representing  no  correction  to  the  model. 

The  and  E  layers  are  separated  by  composition,  i.e.,  E  layer  adjustments  are 

affected  by  adjusting  molecular  ions  (NO^  and  while  layer  adjustments  are  affected 

by  adjusting  the  atomic  ions  (presently  O^,  but  perhaps  including  H"*"  and  He'*’  in  future 
versions).  The  critical  frequencies  if  F^  and  f^E)  are  adjusted  by  scaling  the  appropriate 
ion  density  profiles.  The  layer  heights  adjusted  by  shifting  the 

appropriate  ion  density  profiles  in  altitude.  The  topside  parameters  affect  only  the  (and 

eventually  and  He'*’)  profiles.  The  bottomside  parameters  affect  only  the  molecular  ion 
density  profiles.  The  critical  frequencies,  layer  heights,  and  the  bottomside  parameters  are 
determined  from  DISS  data  and  data  from  the  DMSP  UV  sensors.  The  topside  parameters 
are  determined  from  DMSP  SSIES  and  UV  sensor  data. 

The  eight  adjustment  parameters  are  determined  by  fitting  a  simple  function  of 
latitude  and  longitude  to  the  available  data.  Presently,  PRISM  uses  a  simple  six  parameter 
function  regardless  of  how  much  and  what  kind  of  real  time  data  is  available.  However,  as 
we  gain  experience  using  PRISM  with  a  variety  of  data  sets,  we  will  experiment  with 
additional  parameters  to  determine  the  optimal  number  of  parameters  under  all  possible 
operational  conditions. 

Because  the  real  time  data  is  not  uniformly  distributed  over  the  globe,  there  are 
always  large  regions  devoid  of  data.  The  use  of  SSUSl  images  ameliorates  this  situation 
somewhat,  but  does  noi  completely  eliminate  it.  To  prevent  pathological  values  of  /(A,<p) 
in  data  starved  regions,  we  introduce  "phantom  data"  from  "phantom  sites"  located  in  those 


regions  (Section  4.1).  The  parameter  being  fit  vanishes  at  the.se  phantom  sites  forcing  the 
adjusted  ionosphere  to  "relax"  to  the  model  in  the  data  starved  regions. 

High  Latitudes.  The  previously  developed  high  latitude  specification  model,  HLISM,  has 
been  incorporated  into  PRISM  and  is  used  for  the  trough,  auroral,  and  polar  cap  ionosphere 
with  a  matching  algorithm  to  insure  a  smooth  transition  between  the  midlatitude  and  high 
latitude  regions.  HLISM  uses  a  combination  of  parameterized  physical  models  and 
semi-empirical  models.  Boundaries  between  the  trough,  auroral  oval,  and  polar  cap  regions 
are  determined  from  DMSP  SSIES  and  SSJ/4  data  first.  Then  the  regional  models  are 
adjusted  using  all  available  data.  When  they  become  available,  DMSP  SSUSI  images  will 
be  used  for  both  boundary  determination  and  model  adjustment.  The  details  of  high 
latitude  model  adjustment  vary  depending  on  the  amount  and  kind  of  real  time  data 
available  when  the  model  is  run.  When  data  is  sparse,  the  parameterized  model  is  subject 
to  two  scalings:  first  the  latitude  is  rescaled  according  to  K^,  and  second  a  least  square 
adjustment  is  used  to  provide  an  overall  adjustment  of  When  data  is  abundant,  the 

F  layer  adjustment  makes  use  of  separate  models  for  each  of  the  three  regions  (trough, 
auroral  oval,  and  polar  cap).  The  E  layer  adjustment  uses  measured  characteristics  of 
precipitating  particles  (electrons  and  ions),  a  parameterized  energy  deposition  model,  and  a 
fast  local  chemistry  model  to  calculate  the  E  layer  profile.  Initially,  precipitation  data 
along  the  DMSP  orbital  track  will  come  from  SSJ/4,  but  SSUSI  will  provide  images  from 
which  the  precipitation  parameters  can  be  derived  over  a  large  region. 
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Section  2:  The  Physical  Models 


Four  separate  physical  model  are  used  as  the  basis  of  PRISM:  (1)  a  low  latitude  F 
layer  model  (LOWLAT),  (2)  a  midlatitude  F  layer  model  (MIDLAT),  (3)  a  combined  low 
and  middle  latitude  E  layer  model  (ECSD),  and  (4)  a  high  latitude  E  and  F  layer  model 
(TDIM).  (Since  the  geomagnetic  field  exerts  considerable  control  over  ionospheric 
phenomena,  we  use  the  term  latitude  to  refer  to  magnetic  latitude.)  As  used  to  calculate 
the  PRISM  databases,  all  four  models  use  a  tilted  dipole  magnetic  field  and  corresponding 
magnetic  coordinates.  All  three  models  use  the  MSIS-86  neutral  atmosphere  model  [Hedin, 
1987].  Chemical  reaction  rates,  collision  frequencies,  and  similar  data  are  consistent  among 
all  the  models. 


Section  2.1  The  Low  Latitude  F  Layer  Model 

The  low  latitude  F  region  model  (LOWLAT)  was  originally  developed  by  David  N. 
Anderson  [Anderson,  1973;  reviewed  by  Moffett,  1979).  It  solves  the  diffusion  equation  for 

along  a  magnetic  flux  tube.  Normally,  the  entire  flux  tube  is  calculated  with  chemical 
equilibrium  boundary  conditions  at  both  feet  of  the  flux  tube.  A  large  number  of  flux 
tubes  must  be  calculated  in  order  to  build  up  an  altitude  profile. 

Since  heat  transport  is  not  included,  ion  and  electron  temperature  models  must  be 
used.  We  chose  the  model  of  Brace  and  Theis  [1981].  The  Horizontal  Wind  Model 
(HWM)  of  Hedin  [1988]  was  used  to  describe  thermospheric  winds.  The  critical  feature 
incorporated  in  the  low  latitude  model  is  the  dynamo  electric  field.  The  horizontal 
component  of  this  field  drives  upward  convection  through  the  ExB  drift,  and  this  can 
significantly  modify  profile  shapes  and  densities.  This  phenomenon  is  responsible  for  the 
equatorial  anomaly.  The  ExB  driven  vertical  drift  u.sed  for  these  calculations  was  based  on 
data  from  Jicamarca  [Fejer,  1981,  and  Fejer  et  al.,  1989]  with  modifications  to  simulate 
longitude  variation. 


7 


Section  2.2  The  Midlatitude  F  Layer  Model 


The  midlatitude  F  region  model  (MIDLAT)  is  the  same  as  the  low  latitude  version, 
except  that  the  dynamo  electric  field  is  set  to  zero.  Complete  flux  tubes  are  followed,  but 
neither  horizontal  or  vertical  convection  is  included.  The  computer  resource  requirement  of 
MIDLAT  is  far  less  than  that  of  LOWLAT.  As  long  as  the  boundary  between  low  and 
middle  latitudes  is  chosen  so  that  the  electric  field  is  negligible,  the  two  models  should 
give  identical  results  at  the  boundary  ensuring  continuity  across  the  boundary. 


Section  2.3  The  Low  and  Midlatitude  E  Layer  Model 

The  low  latitude  E  region  model  (ECSD)  was  developed  by  Dwight  T.  Decker  and 
John  R.  Jasperse  and  incorporates  photoelectrons  calculated  using  the  continuous  slowing 
down  (CSD)  approximation  [Jasperse,  1982].  Ion  concentrations  are  calculated  assuming 
local  equilibrium.  A  small  nighttime  source  is  included  to  ensure  a  small  E  layer  is 
maintained  during  the  night. 

Section  2.4  The  High  Latitude  Model 

The  high  latitude  model  (£  and  F  layers)  is  the  Utah  State  University  (USU)  Time 
Dependent  Ionospheric  Model  (TDIM)  developed  by  R.  W.  Schunk  and  J.  J.  Sojka.  [See 
Schunk,  1988  for  a  review.]  This  model  is  similar  to  the  low  and  middle  latitude  models 
except  that  the  flux  tubes  are  truncated  and  a  flux  boundary  condition  is  applied  at  the  top. 
In  addition,  the  flux  tubes  move  under  the  influence  of  the  high  latitude  convection  electric 
field.  In  the  low  latitudes,  because  the  magnetic  field  is  mainly  horizontal,  the  effect  of 
the  electric  field  is  primarily  to  move  the  ionization  in  altitude.  In  contrast,  the  high 
latitude  magnetic  field  is  mainly  vertical,  the  electric  field  driven  convection  is  horizontal. 

The  TDIM  includes  an  f-layer  model  that  incorporates  the  effects  of  ionization  by 
precipitating  auroral  particles.  The  ion  production  rates  used  were  calculated  using  the 
B3C  electron  transport  code  [Strickland,  1976;  Strickland  et  al..  1991}  and  incident  electron 
spectra  representative  of  DMSP  SSJ/4  data.  The  characteristics  of  the  electron  spectra  were 
taken  from  the  Hardy  et  al.  [1987]  electron  precipitation  model. 
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Section  3:  Parameterization  of  the  Physical  Models 

Parameterization  of  the  physical  models  proceeds  in  two  steps.  First,  the  models  are 
used  to  generate  a  number  of  "databases"  for  a  discrete  set  of  geophysical  conditions. 
These  databases  consist  of  ion  density  profiles  on  a  discrete  grid  of  latitudes  and  longitudes 
for  a  24  hour  period  in  UT.  Second,  to  reduce  storage  requirements,  the  databases  are 
approximated  with  "semi-analytic"  functions.  The.se  two  processes  are  described  in  the 
following  subsections. 


Section  3.1  Geophysical  Parameters 

All  the  physical  models  have  been  parameterized  in  terms  of  season  and  solar 
activity.  The  middle  and  high  latitude  models  were  also  parameterized  in  terms  of 
magnetic  activity,  while  the  high  latitude  model  was  additionally  parameterized  in  tenns  of 
the  sign  of  the  interplanetary  magnetic  field  component  (The  high  latitude  model  was 
only  run  using  southward.  Northward  B^  conditions  are  modeled  using  the  low  magnetic 

activity  databases.)  For  the  middle  and  low  latitudes,  the  F  layer  (O"*")  and  the  E  layer 

(NO^  and  O^)  were  computed  and  parameterized  separately.  The  high  latitude  model 
(TDIM)  produced  all  three  ions  simultaneously. 

Due  to  time  and  computer  resource  limitations,  only  a  few  values  of  each  parameter 
are  used.  The  season  "values"  are  the  June  and  December  solstices  and  the  March  equinox 
(which  also  "stands  in"  for  the  September  equinox).  The  values  of  the  other  parameters  are 
summarized  for  each  latitude  region  in  Table  1  on  the  next  page. 


Section  3.2  Representation  of  the  Databases 

When  the  models  are  run  for  any  one  set  of  geophysical  parameters  (e.g.,  June,  ^jq-^ 

=  130.  =  1)-  they  produce  ion  densities  (O^,  NO^,  Op  on  a  grid  of  latitudes,  longitudes, 

and  altitudes  for  a  24  hour  period  of  Universal  Time.  In  order  to  make  this  mass  of 
numbers  more  manageable,  we  produce  a  semi-analytical  representation  of  the  database. 
The  space  and  time  grid  parameters  are  summarized  for  each  latitude  region  in  Table  2. 


9 


Table  1;  Geophysical  Parameter  Values 


Solar  Activity 

F 

10.7 

Magnetic  Activity 

K 

P 

IMF  B 

y 

number  of 

databases 

Low  Latitude 

F  layer 

70,  150,  210 

N/A 

N/A 

36=* 

Midlatitude 

F  layer 

70,  150,  210 

1,  3.5,  6 

N/A 

54b 

Low  &  Midlatitude 

E  layer 

70,  150,  210 

1,  3.5,  6 

N/A 

54c 

High  Latitude 

E  &  F  layer 

70,  130,  210 

1,  3.5,  6 

+,  - 

324d 

33  seasons  X  3  solar  activities  X  4  longitude  sectors 

1*3  seasons  X  3  solar  activities  X  3  magnetic  activities  X  2  hemispheres 

^=3  seasons  X  3  solar  activities  X  3  magnetic  activities  X  2  species 

‘I3  seasons  X  3  solar  activities  X  3  magnetic  activities  X  2  By's 

X  3  species  x  2  hemispheres 


Table  2  Horizontal  Grid  Parameters 


magnetic 

latitude 

magnetic 

longitude 

UT 

number  of 
altitude 
profiles  per 
databa.se 

Low  Latitude 

-32°  to  32° 

0°,  30°.  139°, 

0000  to  2330 

1,584 

F  layer 

in  2°  steps 

and  237° 

in  30  min  steps 

Midlatitude 

30°  to  74°  and 

0°  to  345° 

0100  to  23(K) 

3,456 

F  layer 

-30°  to  -74° 

in  15°  steps 

in  2  hr  steps 

in  4°  steps 

Low  and  Mid- 

-76°  to  76° 

0°  to  345°° 

0100  to  2300 

11,232 

latitude 

in  4°  steps 

in  15°  steps 

in  2  hr  steps 

E  layer 

High  Latitude 

51°  to  89°  and 

7.5°  to  352.5° 

0100  to  2300 

5,760 

E  &  F  layers 

-51°  to  -89° 

in  15°  steps 

in  2  hr  steps 

in  2°  steps 

Due  to  the  computer  resource  requirements  of  the  low  latitude  F  layer  code,  it  was  used  to 
generate  databases  at  four  discrete  longitudes  (corresponding  to  longitude  sectors  for  which 
ExB  drift  measurements  were  available).  Each  longitude  sector  was  parameterized 
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separately,  and  the  necessary  longitude  interpolation  is  carried  out  in  PRISM  during 
execution. 

The  parameterization  of  each  database  was  done  in  three  phases: 

Phase  1.  Empirical  Orthonormal  Functions  (EOF’s)  for  altitude  profiles. 

Phase  2.  Fourier  series  for  magnetic  longitude  variation  (except  low  latitude  F  layer 
where  the  UT  variation  was  fit) 

Phase  3.  Discrete  Orthogonal  Polynomials  for  the  magnetic  latitude  variation 

Except  for  the  low  latitude  F  layer,  all  coefficients  are  tabulated  in  UT  rather  than  fit  with 
analytic  functions.  For  the  low  latitude  F  layer,  the  UT  variation  was  fit  with  Fourier 
series,  but  the  coefficients  are  tabulated  for  the  four  longitude  sectors. 


Section  3.2.1  Phase  1:  Empirical  Orthonormal  Functions 

This  treatment  of  empirical  orthogonal  functions  (EOF.s)  is  based  on  the  Appendix  of 
Secan  and  Tascione  [1984],  which  was  based  on  Lorenz  [1956],  Kutzbach  |1967].  and 
Davis  [1976].  The  reader  is  referred  to  these  references  for  mathematical  proofs  of  the 
assertions  made  below. 

A  database  consists  of  altitude  profiles  at  certain  longitudes,  certain  latitudes,  and 
certain  Universal  Times  (Table  2).  Let  N  =  the  number  of  altitude  profiles  in  a  database, 
and  let  M  be  the  number  of  points  in  each  altitude  profile. 

We  would  like  to  represent  each  altitude  profile  of  the  quantity  yr  as  an  expansion  in 
orthogonal  functions, 

K 

w  (z  )  =  ^  a  ,g,(z  )  +  r  (z  ),  n=l...N, 

m  /  M  nk^k  m  n  m 

k=\ 

where  r  (z  )  is  the  residual  and  the  coefficients  a  ,  are  calculated  from 

n  m  nk 


11 


u 

O  t  ^  ¥  (Z  )  gAz  ) 
nk  nr  m 

m=l 

In  principle,  any  orthogonal  set  of  functions  may  be  used.  However,  the  references  provide 
an  algorithm  for  finding  the  set  which  minimizes  the  RMS  error  for  a  given  number  of 
terms,  K. 

First  define  the  Af  by  Af  covariance  matrix  C  with  elements 


C. 

y 


N 


n=l 


ij  = 


Now  consider  the  eigenvalue/eigenvector  problem  C(f>  =  or 


M 


I 

J  =  l 


M 


I  ¥/. 


^ik  k 


where  ^  is  the  matrix  of  eigenvectors  of  C,  and  L  is  a  diagonal  matrix  whose  elements  are 

th 

the  corresponding  eigenvalues.  (The  k  column  of  ^  is  the  eigenvector  corresponding  to 

th 

the  k  eigenvalue,  By  convention,  the  eigenvectors  and  eigenvalues  are  ordered  so 

that  A|  >  ^2  >  ...  >  Because  C  is  a  real  symmetric  matrix,  eigenvectors  corresponding 

to  unique  eigenvalues  are  guaranteed  to  be  orthogonal  jSee,  e.g.,  Hildebrand,  1965]. 
Because  of  the  origin  of  the  matrix  C,  it  is  unlikely  that  any  of  its  eigenvalues  will  be 
degenerate,  so  we  may  assume  that  ^  is  an  orthogonal  set. 

According  to  Secan  and  Tascione  [1984]  and  references  therein,  the  set  of  orthogonal 
functions  that  minimizes  the  RMS  error  for  K  terms  is  just  the  first  K  eigenvectors; 
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These  are  the  Empirical  Orthogonal  Functions  (EOFs). 

As  a  practical  matter,  we  have  found  that  the  number  of  EOFs  needed  to  provide  a 
reasonably  good  representation  for  all  the  profiles  is  about  Af/6,  as  illustrated  in  Table  3. 

The  only  exception  is  the  low  and  midlatitude  E  layer  (N^  &  Op,  probably  because  these 
databases  covered  both  hemispheres  simultaneously.  We  have  also  found  that  substantial 
improvement  in  representation  does  not  occur  until  the  number  of  EOF's  is  about  A//2. 
Furthermore,  the  EOFs  derived  for  one  database  were  inadequate  for  any  other  database, 
and  the  EOFs  simultaneously  derived  from  several  databases  produce  noticeably  poorer 
representations  than  those  derived  for  each  database  individually.  Consequently,  we  have 
derived  separate  EOF  sets  for  each  database. 

Table  3:  Altitude  Grids  and  EOFs 


Database 

number  of 
alt.  pts. 

minimum 

altitude 

maximum 

altitude 

number  of 
EOFs 

low  latitude  0^ 

55 

160 

1600 

9 

midlatitude  0^ 

49 

125 

1600 

8 

low  &  midlatitude  &  0^ 

28 

90 

400 

7 

high  latitude  O^,  Np  &  Oj 

37 

100 

800 

6 

Note  that  in  none  of  these  cases  was  the  altitude  spacing  unifonn. 


The  six  EOFs  derived  from  the  high  latitude  (USU)  database  for  December,  moderate 
magnetic  activity,  moderate  solar  activity,  and  the  BC  convection  pattern  for  each  of  the 
three  ions  are  shown  in  Figure  1.  The  first  EOF  always  has  the  least  structure,  and 

successive  EOFs  become  more  structured.  The  EOFs  for  low  latitude  and  midlatitude  O"*^ 

profiles  and  for  low  and  midlatitude  NO^  and  O^  are  shown  in  Figure  2. 
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altitude  (km) 


-1.0  0.0  1.0  -1.0  0.0  1.0 


Empirical  Orthonormal  f-unctions  Empirical  Orthonormal  hunctions 


Figure  1:  The  first  six  Empirical  Orthonormal  Functions  for  (a)  0+,  (b)  NO"^,  and  (c) 

O2  derived  from  the  databases  for  December  solstice,  moderate  magnetic 
activity,  moderate  solar  activity,  and  the  BC  (By>0)  Heppner-Maynard 
convection  pattern. 
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altitude  (km)  altitude  (km] 


Empirical  Orthonormal  Functions 


Empirical  Orthonormal  Functions 

Figure  2:  The  Empirical  Orthonoimal  Functions  (EOF’s)  for  (a)  low  latitude  (b) 

midlatitude  O"*^,  (c)  low  and  midlatitude  NO^,  and  (d)  low  and  midlatitude 

derived  from  the  databases  for  December  solstice,  moderate  magnetic 
activity,  and  moderate  solar  activity.  For  (a)  the  USA  longitude  sector  is 


Section  3.2.2  Local  time  and  latitude  representations 


Because  we  were  trying  to  represent  discrete  data  (rather  than  continuous  functions), 
and  because  we  were  working  with  regional  rather  than  global  data  sets,  we  felt  that  the 
standard  spherical  harmonic  expansion  techniques  were  not  appropriate.  Instead,  we  chose 
to  use  orthogonal  functions  of  discrete  variables.  For  longitude  variations  (and  for  the  low 
latitude  F  layer  UT  variation),  the  obvious  choice  was  Fourier  .series,  since  trigonometric 
functions  retain  their  orthogonality  properties  on  unifonn  discrete  grids  and  because  the 
data  is  periodic  in  the  independent  variable.  For  the  latitude  variations,  we  chose  to 
generate  orthogonal  polynomials  using  the  algorithm  derived  in  Beckmann  [1973]  and 
described  below. 

We  proceeded  in  two  stages:  (IJ  fit  the  EOF  coefficients,  with  Fourier  series  at 
each  latitude,  [2]  fit  the  Fourier  coefficients  with  polynomials  in  latitude.  Recall  that  the 
subscript  k  refers  to  the  EOF  while  the  subscript  n  refers  to  a  point  on  the  latitude, 
longitude,  UT  grid.  To  make  this  explicit,  let  us  replace  with  where  A  is 

magnetic  latitude,  (p  is  magnetic  longitude,  and  T  is  Universal  Time.  Let  the  number  of 
latitude  points  be  P,  the  number  of  longitude  points  be  Q,  and  the  number  of  Universal 
Time  points  be  R.  Then  the  total  number  of  points.  A/,  is  PQR  and  the  subscript  ranges  are 
1  <  p  <  P,  0  <  q  ^  Q-\,  and  0  <  r  ^  R-l.  One  way  of  numbering  the  points  is  n  =  p  + 
qP  +  rQP,  but  any  method  that  assigns  a  unique  number  to  each  point  is  acceptable. 

The  Fourier  representation  of  a^(A^,<p^,T^  is 

where 


Q-\ 

q=Q 


a 


kO 


■  -f  * 


J 

1 

7  =  1 
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and 


Q-1 

<7=0 

For  all  of  the  databases  we  found  that  7  =  4  (giving  9  Fourier  coefficients)  reproduced  the 
longitude  behavior  of  the  a^'s  quite  well. 

The  next  step  was  the  representation  of  the  a^ys  and  /)^ys  as  functions  of  magnetic 
latitude.  The  algorithm  for  generating  orthogonal  polynomials  on  a  specified  grid  is  given 
by  Beckmann  [1973].  Let  us  denote  the  desired  polynomials  by  m.(A)  and  define  u  j(A)  =  0 
and  Mq(A)  s  1.  Let  us  denote  the  specified  grid  by  the  set  of  values  (Aj,  A.,,  ...  A^}.  The 
recursion  relation  for  the  polynomials  is 

h^. 

w.^j(A)  =  (A-fipM.(A)  -  ,(A) 

*/-i 

where  the  norms  h.  are  given  by 

P 

/i.  =  y  m^(A) 

p=l 

and  the  recursion  constants  B.  are  given  by 

P 

=  -I  y 

I  ^2  p  I  p 

'  P=1 

The  polynomials  generated  by  this  algorithm  may  be  used  to  represent  the  latitude 
variations  of  the  a, .  and  b, 

kj  kj 
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/ 


VW  -  SW“/‘V 

t=l 

/ 

vw  - 

/=! 

where 

P 

W  '  ^XWV",<V 

'  p=l 

p 

?  Evvv  “,<v 

'  P=1 

The  value  of  I  that  gives  a  good  representation  of  the  latitude  variation  of  the  Fourier 
coefficients  depends  on  the  number  of  latitude  points  in  the  database  grid  and  the  amount 
of  latitude  structure  in  the  Fourier  coefficients. 


low  latitude  F  layer: 

P  =  33 

1  =  12 

midlatitude  F  layer: 

P  =  12 

1  =  9 

low  &  midlatitude  E  layer: 

P  =  39 

I  =  13 

high  latitude  E  &  F  layer: 

P  =  20 

1  =  9 

For  most  of  the  databases,  the  low  order  c,  .  and  d, .  show  a  simple  diurnal  variation 

kji  kji 

in  Universal  Time,  but  the  higher  order  coefficients  show  some  irregular  structure. 
Consequently,  we  have  chosen  not  to  produce  an  analytic  representation  of  the  UT  variation 
and  have  simply  tabulated  the  coefficients  at  the  given  values. 

To  reconstruct  an  O’*"  altitude  profile  in  the  midlatiiude  and  high  latitude  F  layers,  or 
a  molecular  ion  profile  anywhere,  one  must  perform  the  following  steps. 
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(1)  Interpolate  the  orthogonal  polynomial  coefficients  and  in  UT. 

(2)  Reconstruct  the  Fourier  coefficients  for  the  desired  latitude. 

(3)  Reconstruct  the  EOF  coefficients  for  the  desired  longitude. 

(4)  Reconstruct  the  altitude  profile  from  the  EOF’s. 

To  reconstruct  an  O'*’  altitude  profile  in  the  low  latitude  region,  one  must  perform  a 
similar  sequence  of  steps. 

(1)  Reconstruct  the  Fourier  coefficients  for  the  desired  latitude  for  each  of  the  four 
longitude  sectors. 

(?)  Reconstruct  the  EOF  coefficients  for  the  desired  UT  for  each  of  the  longitude 
sectors. 

(3)  Reconstruct  the  altitude  profile  for  each  of  the  longitude  sectors. 

(4)  Interpolate  to  obtain  the  altitude  profile  at  the  desired  longitude. 

The  interpolation  in  step  (4)  is  carried  out  as  a  four  step  process  using  the  following 
interpolation  formula: 

(a)  Use  Fourier  interpolation  on  the  fJF^  and  values  from  each  longitude  sector 
to  obtain  the  values  at  the  desired  longitude. 

(b)  Scale  each  profile  to  match  those  values  (using  the  algorithms  developed  for  the 
real  time  adjustment  process  and  described  in  Section  4.4). 

(c)  Use  Fourier  interpolation  on  the  scaled  profiles  to  obtain  an  interpolated  profile  at 
the  desired  longitude. 

(d)  Rescale  the  interpolated  profile  to  ensure  that  its  and  match  the  values 
obtained  in  step  (a). 

The  reconstructed  profiles  are  tabulations  of  concentrations  on  a  discrete  altitude  grid. 
Interpolation  may  be  u.sed  to  obtain  concentrations  on  altitude  points  lying  between  grid 
points. 
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Section  3.3  Merging  the  Regional  Models 


Because  we  used  four  different  regional  models  in  the  development  of  PRISM,  the 
models  must  be  merged  at  region  boundaries.  In  panicular,  the  low  latitude  and 

midlatitude  O'*’  models  have  to  be  merged  across  the  boundary  between  low  and  middle 

latitudes.  All  three  ions  (O^,  NO^,  O^)  mast  be  merged  across  the  boundary  between 
midlatitudes  and  high  latitudes. 

The  transition  from  low  latitude  profiles  to  midlatitude  profiles  takes  place 
between  30°  and  34°.  The  transition  is  accomplished  by  taking  a  weighted  average  of  the 
h^F^  values  from  the  two  models  in  which  the  weight  shifts  linearly  from  100%  low 
latitude  at  30°  to  100%  midlatitude  at  34°.  The  profiles  are  shifted  to  match  the  averaged 
h^F.^  values  and  then  a  similar  weighted  average  of  the  shifted  profiles  is  taken  to  produce 

the  final  merged  profile.  No  transition  for  NO^  and  is  necessary  since  a  single  model 
was  used  for  these  ions. 

The  transition  from  midlatitude  to  high  latitude  takes  place  over  an  8°  wide  zone 
whose  poleward  boundary  is  the  equatorward  boundary  of  the  trough.  The  transition 
process  is  similar  to  the  low  to  midlatitude  transition,  except  that  the  high  latitude  profiles 
are  shifted  to  match  the  h  F..  and  h  E  values  given  by  the  midlatitude  models.  The  final 
profile  is  produced  by  a  weighted  average  of  midlatitude  and  (shifted)  high  latitude  profiles. 
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Section  4;  Real  Time  Adjustment  Algorithm 


The  Real  Time  Adjustment  (RTA)  algorithm  for  the  low  and  middle  latitude  regions 
is  different  from  the  algorithm  used  in  HLISM.  This  is  partly  due  to  the  relatively 
complex  morphology  of  the  high  latitude  ionosphere  and  partly  due  to  an  evolution  in  our 
ideas  about  the  real  time  adjustment  process  resulting  from  our  experience  in  developing 
HLISM. 

Section  4.1  Available  Data 

The  near  real  time  data  available  for  use  in  the  adjustment  process  comes  from  the 
Digital  Ionospheric  Sounding  System  (DISS),  the  Transionospheric  Sensing  System  (TISS). 
and  a  suite  of  Special  Sensors  on  the  DMSP  satellites.  These  latter  include  the  in  situ 
plasma  properties  measured  by  the  SSIES  instrument,  the  precipitating  panicle 
measurements  of  the  SSJ/4  instrument,  remotely  sensed  electron  density  profiles  (EDP's) 
from  the  SSUSI  (multi spectral  UV  irrager)  and  SSULl  (multispectral  limb  imager) 
instruments.  The  two  UV  imagers  will  be  flown  on  DMSP  beginning  some  time  in  the  last 
half  of  the  decade. 

The  DISS  network  consists  of  ahou,  20  Digisondes  located  mainly  in  the  northern 
hemisphere.  These  digital  iono.sondes  return  parameters  that 

allow  the  reconstruction  o^  a  bottomside  profile.  The  DISS  sites  (as  of  15  June  1990)  are 
shown  in  Table  4. 

The  TISS  network  consists  of  about  17  suitably  equipped  GPS  receivers  located 
mainly  in  the  northern  hemisphere,  a  number  of  which  are  co-located  or  nearly  co-located 
with  DISS  sites.  The  projected  TISS  sites  (as  of  9  November  1989)  are  shown  in  Table  5. 
These  receivers  will  provide  Total  Electron  Content  (TEC)  between  the  site  and  the  GPS 
satellites  (about  1 1  earth  radii  in  altitude).  When  the  full  constellation  of  GPS  satellites  is 
operational,  there  will  always  be  at  lea.st  four  .satellites  in  view  and  .sometimes  as  many  as 
eight.  Thus,  considerable  information  about  the  ionosphere  will  be  produced. 
Unfortunately,  because  the  quantity  provided  is  an  integral  quantity  (TEC),  it  is  not  easy  to 
incorporate  it  into  the  PRISM  real  time  scheme.  Therefore,  we  have  decided  not  to  include 
TISS  data  in  Version  1.0.  However,  we  have  given  development  of  a  TISS  algorithm  first 
priority  for  code  improvements  to  be  undertaken  after  Version  1.0  has  been  validated. 
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Table  4:  DISS  Sites 
(15  June  1990) 


site  Name  ICAO  WMO#  lat  Ion 


a 

Amchitka,  AK 

PAHT 

70454 

51.38 

179.37 

b 

Argentia,  NF 

CYAR 

71807 

47.25 

-54.00 

c 

Bermuda 

TXKF 

78016 

32.37 

-64.68 

d 

Bradshaw  AAF,  HI 

PHSF 

91197 

19.78 

-153.55 

e 

College,  AK 

PAEI 

70265 

64.91 

-147.93 

f 

DiyaAakir,  TU 

LTCC 

17280 

37.88 

40.18 

g 

Dyess  AFB,  TX 

KDYS 

72266 

32.40 

-99.86 

h 

Eglin  AFB,  FL 

KVPS 

72221 

30.50 

-86.50 

i 

Goose  Bay,  LB 

CYYR 

71816 

53.32 

-60.43 

j 

Learmonth,  AS 

APLM 

94302 

-22.33 

114.03 

k 

Manila,  RP 

RPMM 

98429 

14.70 

121.10 

1 

McClellan  AFB,  CA 

KMCC 

72483 

38.67 

-121.40 

m 

Petersburg,  AK 

PPSG 

70386 

56.80 

-132.92 

n 

Osan  AB,  Korea 

RKSO 

47122 

37.10 

127.03 

0 

Ramey,  PR 

TJFF 

78514 

18.50 

-67.17 

P 

Sondrestrom,  GL 

BGSF 

04231 

68.00 

-51.00 

q 

Trieste,  Italy 

LIVT 

16110 

45.65 

13.75 

r 

Vandenberg  AFB,  CA 

KVBG 

72393 

34.76 

-120.56 

s 

Wallops  Is,  VA 

KWLI 

72402 

37.93 

-75.47 

t 

RAF  Wethrsfield 

EGVT 

03688 

51.97 

0.50 

site 

Table  5:  TISS  sites 
(9  November  1989) 

Name  WMO#  lat 

Ion 

alt 

(ft) 

a 

Otis  ANGB,  MA 

72506 

41.67 

-70.50 

130 

b 

Shemya  AFB,  AK 

70414 

52.67 

174.17 

100 

c 

Oahu,  HI 

91182 

21.50 

-158.00 

840 

d 

Goodfellow  AFB,  TX 

72263 

31.33 

-100.50 

1895 

e 

Goose  Bay,  CN 

71816 

53.33 

-60.50 

155 

f 

College,  AK 

70265 

64.83 

-147.83 

500 

g 

RAF  Croughton,  UK 

03655 

52.00 

-1.25 

440 

h 

Thule  AB,  GL 

04205 

76.50 

-68.75 

250 

i 

Diego  Garcia 

61967 

-7.33 

72.42 

10 

j 

Howard  AB,  PN 

78806 

8.92 

-79.58 

50 

k 

Guam 

91217 

13.58 

144.92 

635 

1 

Churchill,  CN 

71913 

57.50 

-94.17 

95 

m 

Keflavik,  Iceland 

04018 

64.00 

-22.50 

170 

n 

Ramey,  PR 

78514 

18.50 

-67.17 

250 

0 

Beale  AFB,  CA 

72483 

39.13 

-121.42 

115 

P 

Diyarbakir,  TU 

17280 

37.83 

40.17 

2250 

q 

Kunsan  AB,  ROK 

47140 

35.92 

126.67 

50 
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The  DMSP  satellites  are  in  approximately  circular  orbits  near  840  km  altitude  with 
inclinations  chosen  to  keep  the  orbit  "sun  synchronous”  -  i.e.,  the  local  time  of  the  orbit  is 
approximately  constant.  There  are  usually  two  operational  satellites  at  any  time,  one  with 
an  orbit  near  0430-1630  local  time  and  the  other  near  0730-1930  local  time.  Current 
satellites  carry,  among  other  instruments,  the  SSIES  instrument,  which  returns  ion  and 
electron  densities  and  temperatures,  ion  drift  velocities,  and  the  SSJ/4  instrument,  which 
measures  precipitating  differential  ion  and  electron  fluxes  between  20  eV  and  20  keV. 
Beginning  near  the  end  of  the  decade,  DMSP  satellites  will  also  carry  two  UV  imagers, 
SSUSI  looking  down  and  SSULI  looking  at  the  limb.  The  airglow  and  auroral  optical 
intensities  measured  by  these  two  instruments  will  be  used  to  determine  the  electron  and 
neutral  density  profiles.  It  is  not  yet  clear  whether  the  EDP's  will  be  provided  to  PRISM, 
or  whether  the  eight  parameters  used  in  PRISM's  profile  adjustment  algorithm  will  be 
provided.  Either  way,  PRISM  will  ultimately  use  tho.se  eight  parameters  rather  than  the 
full  EDP's.  Initially,  each  UV  instrument  will  provide  EDP's  (or  EDP  parameters)  based  on 
its  own  observations  only  (single  instrument  EDP's).  However,  it  is  hoped  that  eventually 
EDP's  will  be  derived  from  all  four  instruments  (SSUSI,  SSULI,  SSIES,  and  SSJ/4).  These 
should  be  much  more  accurate  because  the  data  from  these  instruments  is  complementary. 


Section  4.2  Low  and  Midlatitude  Adjustment  Parameters 

The  RTA  will  operate  on  eight  parameters  that  prescribe  how  an  electron  density 
profile  is  to  be  modified  or  "corrected"; 

(1)  AfF.,,  the  correction  to  the  model /F.,  calculated  from  O''^  profile 

(2)  correction  to  the  model  /^£,  calculated  from  NO^  and  profiles 

(3)  A/?  F-,  the  correction  to  the  model  h  F-,  calculated  from  the  O'*’  profile 

m2  m2  ^ 

(4)  Ah  E,  the  correction  to  the  model  h  E,  calculated  from  the  NO'*’  and  Ot 

m  m  2 

profiles 
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(5)  ^lop’  correction  to  the  O"*^  density  at  a  specific  altitude  (i.e.,  the  DMSP 
altitude) 

(6)  ^top'  correction  to  the  scale  height  at  a  specific  altitude  (i.e.,  the 
DMSP  altitude) 

(7)  A/Vj,  the  correction  to  the  electron  density  8  km  below  hJE,  and 

(8)  A/Vj,  the  correction  to  the  electron  density  16  km  below  h^E. 

The  nominal  value  for  each  of  these  parameters  is  zero.  A  positive  (negative)  value  means 
that  the  model  value  must  be  increased  (decreased).  Using  the  available  near  real  time 
data,  the  RTA  process  will  assign  non-zero  values  as  functions  of  magnetic  latitude  and 
longitude.  The  RTA  process  is  described  below.  The  use  of  the  eight  parameters  to 
modify  the  model  ion  density  profiles  is  described  in  Section  4.4. 


Section  4.3  Real  Time  Adjustment  of  the  Low  and  Midlatitude  Profile  Parameters 

In  PRISM  Version  1.0,  Parameters  (1),  (3),  (5),  and  (7)  apply  to  the  profile  only, 

while  the  remaining  four  apply  to  the  molecular  ion  (NO"'’  and  O^)  profiles  only.  In  either 
case,  the  critical  frequency  and  layer  height  adjustment  occurs  first,  followed  by  the  topside 
or  bottom  side  adjustment,  as  appropriate.  Each  of  the  eight  adjustment  parameters  are 
determined  by  fitting  a  simple  function  of  latitude  and  longitude  to  the  available  data.  The 
nominal  functional  form  is 

y(A,<p)  =  Pj  +  +  PjCos(2A)  +  /7^sin(3A) 

+  +  p^sinA  +  p.^cos(2A)  +  Pj^sin(3A)l  simp 

+  IPq  +  P,Q.sinA  +  PjjCos(2A)  +  P|.,sin(3A)|  coscp 

with  twelve  parameters.  When  used  with  data  from  the  20  station  DISS  network  only,  we 
found  it  best  to  set  p^,  Py  Pg,  Pj^,  Pjj,  and  Pp  to  zero,  resulting  in  a  simpler,  six 
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parameter  function. 


/(A,(p)  =  Pj  +  p^sinX  +  p^cosilX)  +  p^sin(3A)  +  p^sincp  +  PgCos(p 

When  SSUSI  data  becomes  available,  potentially  providing  thousands  of  EDP's  rather  than 
20,  the  additional  parameters  (or  perhaps  even  more)  can  be  used.  However,  for  simplicity 
and  consistency,  we  have  retained  the  six  parameter  representation  even  when  using 
simulated  SSUSI  data.  As  we  gain  experience  u.sing  PRISM  with  a  variety  of  .simulated 
data  sets,  we  will  determine  the  optimum  parameter  set  to  use  with  a  given  complement  of 
data.  Due  to  the  modular  design  of  PRISM,  it  will  be  easy  to  upgrade  the  function  fitting 
algorithm. 

Because  the  real  time  data  is  not  uniformly  distributed  over  the  globe,  there  are 
always  large  regions  devoid  of  data.  This  is  especially  true  when  SSUSI  data  is  absent. 
To  prevent  pathological  values  of  f(X,(p)  in  data  starved  regions,  we  introduce  "phantom 
data"  from  "phantom  DISS  sites"  located  in  those  regions.  The  parameter  being  fit 
vanishes  at  these  phantom  sites  forcing  the  adjusted  ionosphere  to  "relax"  to  the  model  in 
the  data  starved  regions.  (Because  the  fitting  process  is  a  least  squares  adjustment,  the 
adjustment  do  not  relax  exactly  to  zero,  but  rather  to  small  values.)  The  phantom  sites  are 
"hard  wired"  at  the  locations  given  in  Table  6,  which  are  only  used  when  SSUSI  data  is 
absent.  So  far,  we  have  not  found  it  necessary  to  use  phantom  data  when  SSUSI  data  is 
available  because  the  extensive  spatial  coverage  provided  by  SSUSI  images  does  not  leave 
large  areas  uncovered.  However,  we  have  assumed  that  SSUSI  will  provide  all  eight 
profile  parameters  for  each  complete  orbit.  It  remains  to  be  determined  whether  or  not  the 
SSUSI/SSULI  topside  parameters  for  the  daytime  ionosphere  will  be  sufficiently  accurate 
for  use  in  PRISM. 


Section  4.4  Modifying  the  Low  and  Midlatitude  Model  Profiles 

The  eight  parameters  will  be  used  to  adjust  the  ion  densities  profiles  in  the  following 

ways. 

Parameters  (1)  and  (2),  the  critical  frequencies,  will  be  used  to  scale  the  entire 
profile.  The  F  layer  (O^)  algorithm  is 
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max 
0<  z<« 


[[0^](z) 


where  [O'^’Kz)  is  the  O'*’  concentration  (cm'^)  at  altitude 
z. 


(b) 


/^F2(model)  =  8.98  x  10'^ 


N  F.(model),  N  F.  in  cm’^,  /  F.  in  MHz 
ml  ml  o  1 


(c)  /^F^Ccoirected)  =  /^F^Cmodel)  + 

(d)  N  F-(corrected)  =  1.24  x  10“^  (f  F.)^ 

ml  0  1 


N  F,  (cor  rected) 

(e)  [0'''](z, corrected)  =  — - (0^](z,model), 

//F(  model) 

m  z 


and  the  E  layer  algorithm  is 


(a) 


E 

m 


max 
0<  z<«> 


[NQ-^Kz)  +  lO^I(z) 


where  jNO^Kz)  and  [O^Kz)  are  the  NO^  and 
©2  concentrations  (cm’^)  at  altitude  z 


(b) 


/^£(model)  =  8.98  x  10*^ 


N  £(model). 
m 


(c)  /^£(corrected)  =  /^£(modeI)  +  /V^£ 


(d)  N  £(corrected)  =  1.24  x  10'‘  (f  E)^ 

m  o 
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N  E (cot  rected) 

(e)  [NO  ](z,corrected)  =  — - [NO'^l(z,modeI), 

N^£(  model) 

N  £( cor  rected) 

[0_](z,corrected)  =  — - [0*](z,model), 

^  N  £(model)  ^ 

m 


Parameters  (3)  and  (4),  the  layer  heights,  will  be  used  to  shift  the  ion  density  profiles 
in  altitude.  These  corrections  operate  on  profiles  that  have  already  been  corrected 
according  to  parameters  (1)  and  (2). 

[0''^l(z,corrected)  =  [0^](z-A/i  £  .model) 

m  2. 

[NO'‘’](z,corrected)  =  [NO'‘’J(z-A/i  £,model) 

tn 


[0!tl(z,corrected)  =  [Ot](z-Ah  £, model) 

X  /.  rn 


Parameters  (5)  and  (6)  are  used  to  correct  the  topside  profile  based  on  n^,  n.,  T^,  and 
T.  measurements  from  SSIES  on  DMSP  (nominally  840  km).  Parameters  (7)  and  (S)  are 
used  to  correct  the  bottomside  profiles  based  on  bottomside  profiles  from  the  DISS 
network.  These  four  parameters  are  determined  (0er  parameters  (1)  through  (4)  have  been 
determined,  and  these  corrections  are  applied  to  profiles  that  have  already  been  corrected 
according  to  parameters  (1)  through  (4).  The  topside  correction  using  parameters  (5)  and 
(6)  follows. 


Let  N  (z)  be  the  model  profile  (after  f  £.  and  h  £»  corrections  have  been 
fn  ^  ■'  o  2  m  2 

applied).  Let  N  (z)  be  the  corrected  profile  based  on  AN  and  A//  Further,  let 
c  ^  lop  lop 

1  =  A  £-  and  z  .  =  the  altitude  at  which  AN,  and  A//  were  determined.  Then  define 
p  rn  l  a  top  lop 


X 


z 


z 

_p_ 

z 

p 
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so  that 


The  bottomside  correction  proceeds  as  follows.  Let  be  the  model  NO"^  + 

profile  {after  fJE  and  h^E  corrections  have  been  applied).  Let  A/ (2)  be  the  corrected 
profile  based  on  bottomside  data.  Further,  let  and  let  Z|  and  be  two  altitudes 

at  which  the  actual  electron  density  has  been  measured.  Then  the  corrections  to  at 

those  altitudes  are 
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ANj  =  -  N^(2) 

For  convenience,  choose  Zj  =  (z^  +  22)^-  Then  define 

JC  = - 

‘  ^0 

(Note  that  jr  is  negative  for  -«  <  z  <  z^.) 

>>^U)  =  In  N^(z) 

y jx)  =  In  N  (z) 
m  w 

g{x)  =  -  yj.x) 

so  that 

A/^(z)  =  NJz)  explgCr)] 

Further  define 

/?j  s  In  1  + 
and 

/?2  =  In  1  + 

Let  g(x)  have  the  form 

0. 

gj(jc),  -2  <  or  <  0 

g|(-2),  -~  <  jr  <  -2 

with  gj(x)  =  (/?,  -  ^2^^"*  +  ■  ^2^^^ 


A/Vj  ■ 

N  (z.) 
m  1  ■* 

AA/^  ■ 

^W- 
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This  form  ensures  that  the  correction  and  its  first  derivative  vanish  at  the  peak  {x  =  0),  that 
the  correction  has  the  desired  values  at  z  =  Zj  and  z  =  z^,  and  that  the  correction  is  well 
behaved  as  z  -* 


Section  4.5  The  High  Latitude  Adjustment  Algorithm 

Due  to  the  complexity  of  the  high  latitude  ionosphere,  the  real  time  adjustment 
algorithm  differs  appreciably  from  the  low  and  midlatitude  algorithm.  Until  SSUSI  data 
becomes  available,  there  will  be  insufficient  data  to  adjust  the  parameterized  USU  model  in 
the  way  that  the  parameterized  low  and  midlatitude  models  can  be  adjusted.  Even  with 
SSUSI  data,  it  is  not  clear  that  the  midlatitude  algorithm  is  appropriate  for  the  high  latitude 
regions. 

In  PRISM  Version  1.0,  the  first  step  in  the  high  latitude  real  time  adjustment  process 
is  the  establishment  of  boundary  locations.  Three  boundaries  are  required:  (1)  the 
equatorward  edge  of  the  trough,  (2)  the  equatorward  edge  of  the  auroral  oval,  and  (3)  the 
equatorward  edge  of  the  polar  cap. 

The  equatorward  edge  of  the  trough  is  determined  from  SSIES  drift  meter  data.  The 
trough  boundary  as  a  function  of  magnetic  longitude  is  given  by  the  formula 


0^i(p)  =  0Q(/p  +  a  exp 


where  (p  is  the  magnetic  local  time  (MLT,  hours),  and  is  a  "trough  index"  correlated  with 
K^.  is  the  radius  of  the  trough  boundary  at  magnetic  local  midnight  and  is  given  by 


=  24.4°  +  2.12°/^ 


The  second  term  represents  the  dayside  distortion  of  the  boundary,  which  would  otherwise 
be  a  circle  centered  on  the  magnetic  pole.  The  parameter  values  are 


a  =  -10.5° 
b  =  11.5  hr 
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c  =  3.88  hr 
p  =  2.73 

This  model  is  based  on  the  convection  boundaries  shown  in  Heppner  and  Maynard  [1987], 

The  equatorward  edge  of  the  auroral  oval  is  determined  from  an  analytic 
representation  of  the  Gussenhoven  et  al.  (1983}  boundary.  The  boundary  itself  is  a  circle 
whose  center  is  displaced  from  the  magnetic  pole.  The  radius  (in  degrees)  of  the  circle  is 

=  20.9°  +  1.7°  Pj 

where  PI  is  a  "precipitation  index".  The  center  of  the  circle  is  located  at 

A.  ,(P,)  =  87.3°  -  0.267°  P, 

cl  1  1 

=  39.5°  -  1.25°  Pj  +  0.076°  p\ 

where  is  magnetic  atitude  and  is  magnetic  local  time  (in  degrees). 

The  equatorward  edge  of  the  polar  cap  l  as  almost  the  same  form,  except  that  it  is 
parameterized  in  terms  of  a  separate  "precipitation  index"  P^. 

e^iP^)  =  13.4°  +  1.7°  P^ 

■  89.2°  +  0.267  P^.  P2  <  3 

90.8°  -  0.267°  P^,  ^2  -  ^ 
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Note  that  the  form  of  <p^2  guarantees  that  the  centers  of  the  two  circles  lie  on  a  great  circle 
passing  through  the  magnetic  pole. 

The  next  step  in  the  high  latitude  adjustment  process  depends  on  the  amount  and 
kind  of  data  available.  The  decision  of  how  to  proceed  is  made  separately  for  the  E  layer 

(NO'*’  and  Op  and  the  F  layer  (O'*’).  In  each  case,  two  choices  are  available: 

F  layer 

1.  Perform  a  simple  least  squares  adjustment  of  the  USU  model. 

2.  Use  a  semi -empirical  model  (FMODEL)  to  adjust  the  USU  O'*’  profiles. 

E  layer: 

1.  Perform  a  simple  least  squares  adjustment  of  the  USU  NO'*'  and  models. 

2.  Use  a  fast,  first  principles,  E  layer  local  chemistry  model  (HLE). 

The  decision  is  made  in  PRISM  as  follows: 


DISS 

data  source 
SSIES 

SSJ/4 

model  used 

F  layer  E  layer 

Yes 

Yes 

Yes 

FMODEL 

HLE 

Yes 

Yes 

No 

FMODEL 

USU 

Yes 

No 

Yes 

FMODEL 

HLE 

Yes 

No 

No 

USU 

USU 

No 

Yes 

Yes 

USU 

USU 

No 

Yes 

No 

USU 

USU 

No 

No 

Yes 

USU 

USU 

No 

No 

No 

USU 

USU 

In  experimenting  with  high  latitude  data,  we  found  that  extrapolating  SSJ/4  data  taken 
along  the  DMSP  orbital  track  to  points  well  away  from  the  orbital  track  was  very  risky. 
We  found  no  suitable  model  of  the  instantaneous  auroral  precipitation  for  this  extrapolation. 
When  SSUSI  auroral  image  data  becomes  available,  this  limitation  will  be  removed  because 
much  less  extrapolation  will  be  required.  It  should  be  possible  to  use  HLE  whenever 
timely  SSUSI  images  are  available. 


FMODEL  is  a  semi -empirical  model  of  ^  combination  of  theory  and 

data.  It  is  divided  into  three  regions:  the  subauroral  trough,  the  auroral  oval,  and  the  polar 
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cap.  The  determined  by  least  squares  adjustment  of  the  model  parameters  is  used  to 
scale  the  USU  O'*’  profiles.  No  further  adjustment  of  the  profiles  is  performed. 

The  subauroral  trough  is  divided  into  two  local  time  regimes;  evening  (from  1200  to 
0000  MLT)  and  morning  (from  0000  to  1200  MLT).  Each  trough  (morning  and  evening) 
has  a  depth  parameter  that  specifies  the  difference  between  the  midlatitude  value  of  foF2  at 
the  equatorward  edge  and  the  trough  minimum.  The  local  time  variation  of  the  depth  is 
fixed  (not  part  of  the  least  squares  adjustment  process).  If  the  width  of  the  trough  is  less 
than  3°,  the  depth  is  reduced  in  proportion  to  the  width  so  that  when  the  width  vanishes  so 
does  the  depth.  The  thickness  of  the  (equatorward)  trough  wall  is  always  60%  of  the  total 
width  of  the  trough.  The  poleward  edge  of  the  trough  is  the  equatorward  edge  of  the 
auroral  F-layer,  so  the  poleward  "wall"  is  considered  to  be  part  of  the  auroral  region. 

The  auroral  F  layer  is  simply  a  cubic  polynomial: 

/  =  f  +  A  {I  -  X  f  +  B  (X  -  k  )^ 

•'o  2  ■'max  max  max 

where  /  -  f  FAX  )  is  an  extremum,  and  A  and  B  are  chosen  so  that  /  F,  is 

■’max  'o  2  max  o  2 

continuous  across  the  boundaries  with  the  trough  and  polar  cap: 

The  background  polar  cap  is  obtained  from  the  URSI  coefficients  using  an 
effective  sunspot  number.  This  is  not  the  same  SSNeff  used  with  ICED.  Instead,  the 
value  of  the  polar  cap  effective  sunspot  number  is  determined  as  part  of  ULISM's  least 
squares  adjustment  process  based  on  ionosonde  data.  It  should  be  noted  that  this  part  of 
the  model  describes  only  the  background  polar  cap  ionosphere. 
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Section  S:  PRISM  Execution 


PRISM  has  been  delivered  to  AWS  on  a  DEC  TK50  tape  cartridge  in  VAX  VMS 
BACKUP  format.  Instructions  for  restoring  the  files  accompanied  the  tape.  The  files  on 
the  tape  cartridge  include  the  complete  source  code  for  PRISM,  command  files  for 
compiling,  linking,  and  executing  PRISM,  the  necessary  input  and  data  files  for  execution, 
and  the  output  files  that  resulted  when  PRISM  was  executed  on  CPI/Boston's  MicroVAX. 

The  command  file  [PRISM. COM_FILES]OPTFOR. COM  may  be  submitted  to  fully 
run  the  code  and  reproduce  the  output  file  that  is  included  on  the  tape.  This  command  file 
compiles  all  the  modules  and  calls  OPTLNK.COM.  OPTLNK  links  the  modules  to  produce 
the  executable  code.  OPTLNK  then  calls  the  command  file  88113S_UF,  which  runs  the 
case.  As  supplied  on  the  tape,  this  command  file  requires  that  PRISM  be  the  main 
directory,  and  that  the  directory  structure  be  as  follows. 


If  PRISM  cannot  be  the  main  directory,  or  the  directory  structure  cannot  be 
preserved,  then  the  command  files  must  be  modified  as  follows: 

If  only  SIM_DATA,  GCXCGM,  SGI,  or  DATABASES  is  to  be  changed,  then  edit 
the  file  PATH_NAMES.TXT  and  give  the  new  directories  where  the  respective  data  will 
be. 
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If  the  FORTRAN  modules  of  PRISM  reside  in  a  different  directory,  then  change 
OPTFOR,  OPTLNK  and  88113s_UF  to  give  the  correct  location  of  the  FORTRAN  code. 

If  PRISM  cannot  be  the  main  directory,  all  three  command  files  (OPTFOR,  OPTLNK 
and  88113S_UF  must  be  changed  to  show  the  new  path  to  PRISM.  PATH_NAMES.TXT 
must  also  be  changed. 


Section  5.1  Code  Structure 


A  complete  structure  chart  for  PRISM  appears  on  the  following  six  pages.  Because 
of  the  complexity  of  the  code,  it  is  broken  into  major  modules.  The  top  level  modules  are; 


GETDAT: 


INIT; 

USER_INPUT: 

MIDLAT: 

HLISM: 

OUTPUT: 


which  obtains  the  data  used  to  drive  the  model,  determines  boundary 
locations,  and  makes  certain  algorithm  choices  based  on  the  kind  and 
quantity  of  data  available, 
which  initializes  various  algorithms. 

which  prompts  the  user  for  certain  choices  concerning  code  operations  and 
reads  in  the  necessary  model  databases. 

which  handles  the  real  time  adjustment  of  the  low  and  midlatitude  models, 
which  handles  the  real  time  adjustment  of  the  high  latitude  models, 
which  writes  the  results  (EDP’s)  of  the  adjustment  to  a  file  for  subsequent 
display. 


In  addition,  there  are  two  major  modules  that  service  other  modules: 


REGMOD:  which  is  called  by  MIDLAT,  HLISM,  and  OUTPUT  to  produce  ion  density 

profiles,  scaled  and  adjusted  as  necessary,  using  whichever  model  is 
appropriate  for  the  latitude  and  longitude  in  question. 

HLEMODEL:  which  is  called  by  REGMOD  to  calculate  the  high  latitude  E  layer  from 
first  principles  when  sufficient  data  is  available. 
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PRISM 
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Section  5.2  Input  and  Output 


All  input  files  necessary  for  the  execution  of  PRISM  have  been  supplied  with  it.  No 
additional  user  input  is  required  if  the  supplied  command  file  is  used  to  execute  the 
program. 

Output  consists  of  a  single  large  file  giving  electron  density  profiles  (EDP’s)  on  a 
grid  of  2°  increments  in  latitude  and  5°  increments  in  longitude.  The  grid  covers  all 
latitudes  (-90°  to  90°)  and  longitudes  (0°  to  360°). 


Section  5.3  Sample  Re.sult.s  Using  Simulated  Data 

We  pre.sent  one  .sample  run  of  PRISM.  A  complete  set  of  DISS,  SSIES,  and 
combined  SSUSI  and  SSULI  data  for  the  midlatitude  region  was  simulated  using  FAIM 
(the  Fully  Analytic  Ionospheric  Model).  The  sunspot  number  was  set  to  70  (^jq^  =  120), 
and  was  set  to  5.  High  latitude  data  was  actual  data  used  in  validating  HLISM.  We 
present  here  only  the  midlatitude  results. 

In  Figures  3  and  4,  we  show  contours  of  FAIM.  In  Figures  5  and 

6  we  show  contours  of  fJF^  and  f  E  produced  by  PRISM  based  on  the  simulated  data. 
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PRISM  FOE  Year=1988 


F10P7=  120.0 

KP-  5.00 

SSN  -  70.00 


Contours  at  [0.38,  1.18,...,  7.60] 


Figure  6  Contours  of  f^E  produced  by  PRISM  using  input  data  simulated  by  FAIM 

for  the  date  and  conditions  of  Figure  4 
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